Anomalous relaxations and chemical trends at III-V nitride non-polar surfaces 
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Relaxations at nonpolar surfaces of III-V compounds result from a competition between de- 
hybridization and charge transfer. First-principles calculations for the (110) and (1010) faces of 
zincblende and wurtzite A1N, GaN and InN reveal an anomalous behavior as compared with ordi- 
nary III-V semiconductors. Additional calculations for GaAs and ZnO suggest close analogies with 
the latter. We interpret our results in terms of the larger ionicity (charge asymmetry) and bonding 
strength (cohesive energy) in the nitrides with respect to other III-V compounds, both essentially 
due to the strong valence potential and absence of p core states in the lighter anion. The same 
interpretation applies to Zn II- VI compounds. 
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I. INTRODUCTION 

The III-V nitrides GaN, A1N, and InN are of enormous 
current interesto in blue optoelectronics and high-power 
devices technology. Among the relevant problems in this 
area, there is the high density of threading dislocations 
and domain boundary defects occurring during growth. 
These boundaries often coincide geometrically with the 
non-polar surfaces of the material, so that accurate char- 
acterizations of these surfaces are of interest, and first- 
principles calculations in this area are timely. Although 
wurtzite nitrides are usually grown along the (0001) po- 
lar direction, other possible growth orientations are being 
examined, such as the nonpolar (1010) and (1120) sur- 
faces. Also, thin films of zincblende GaN were grown 
on various substratesa typically along (110) (one of the 
cleavage faces of zincblende). 

Earlier works suggest that the nitrides behave quite 
differently than the ordinary III-V semiconductors such 
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as GaAs or GaP in sever, 
gap-cohesive energy relatian 

dielectrica and piezoelectrica constants. Recent work: 
have pointed out the unusual surface relaxations of GaN 
as a further point of difference. The latter "anomaly" 
would reflect a stronger ionic character of GaN, making 
it similar to the II- VI oxide ZnO, commonlyi-considered 
as highly ionic. First-principles calculationsHLO for the 
(1010) surface of ZnO gave smaller rotations and larger 
contractions than in GaAs and in other II-VI's (in view 
of the similar morphology and electronic structure of the 
(110) and (1010) surfaces, considerations about the relax- 
ation mechanism are quite valid for both). In the present 
paper we take up this problem for the nitrides, studying 
the wurtzite (1010) and zincblende (110) surfaces of GaN, 
A1N, and InN, and reexamining the properties of the ho- 
mologous surfaces of ZnO and GaAs as reference systems. 
We discuss the results in terms of increased ionicity and 
increased cohesive energy of the nitrides as compared to 



other III-V compounds. Our interpretation also fits the 
situation of II- VI compounds, and is compatible with the 
presence of the "anomaly" only for O and N compounds. 

The present first-principJes calculations are based on 
density functional theorytll in the local density approxi- 
mation (LDA) for the exchange-correlation energjj func- 
tional, for which we adopt the Ceperley-Aidenla form 
as parametrizecL— by Perdew and Zunger.E3 Ultrasoft 
pseudopotentialsE-3 have been employed for all the ele- 
ments involved in the calculations. A plane- wave basis 
is used with a cutoff at 25 Ry. For Ga, In, and Zn, wc 
explicitly include the semicore d electrons in the valence. 
Slab supercells were used to simulate the surfaces. The 
results presented here for zincblende (wurtzite) surfaces 
were obtained with symmetric slabs encompassing 8 (9) 
layers, i.e. 16 (18) atoms, whereby all atomic coordinates 
were relaxed to obtain forces below 1 mRy/bohr. A mesh 
of 10 irreducible special k-points (obtained by downfold- 
ing the bulk mesh) is used for both the zincblende and the 
wurtzite surface. All calculations are performed at the 
theoretical lattice constants: a=6.00 bohr, c/a=1.613, 
u=0.376 for ZnO; a=10.60 bohr for GaAs; a=5.81 bohr, 
c/a=1.619, u=0.380 for A1N; a=6.04 bohr, c/a=1.634, 
u=0.376 for GaN; a=6.66 bohr, c/a=1.627, u=0.377 for 
InN (see Ref. ||for details on the optimization procedure). 

The results for the structural parameters of the 
zincblende and wurtzite surfaces are presented in Sec. 
H and Sec. Ill respectively; those for the energetics and 
electronic structure are discussed in Sec. [fv|. In Sec. |v| 



we discuss our findings on the basis of simple chemical 
concepts. 



II. (110) ZINCBLENDE SURFACES 

The relaxations typical of the (110) surface of 
most Ill-V-aiid II- VI compounds have been generally 
interpretedllaO as driven by a loss of sp 3 hybridization 
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towards anion p and cation sp 2 — like character. Upon 
cleavage, charge is transfered from the cation dangling 
bond into the anion dangling bond. The plane contain- 
ing each anion-cation chain running along [110] rotates 
with respect to the ideal surface (see the sketches in Fig.[l] 
and Fig. ||). In each surface dimer, the cation shifts 
downwards, so as to lay nearly in the plane of its three 
neighbors, and rehybridizes to sp 2 -like. The anion shifts 
upward and is bound to its neighbors by p-like back- 
bonds, while it fills up its low-laying s-like state. 

The relaxations are usually expressed via a combina- 
tion of the layer rotation angle 8, the bond rotation angle 
to, and the bond contraction Cb (see Fig.||). Notice that 

9 and u> are independent parameters, since the dimers can 
stretch or shorten besides rotating. Only if Cb=0, 9 and 

10 are related by -s/3 sin lo = sin 6 . In Table | we list our re- 
sults for A1N, GaN, InN, and GaAs, in comparison with, 
the data of Ref. R and (only for GaAs) experiments.E£l 
The values for GaAs are in very good agreement with ex- 
periment. For nitrides, we do indeed confirm an anoma- 
lous behavior: the rotation angles are nearly a half than 
for GaAs, and the bond contractions are appreciable, as 
opposed to negligible for GaAs. If (see the discussion 
in Sec. |v|) we interpret small bond rotations and large 
bond contractions as a measure of ionicity, we see that 
the latter grows along the sequence InN— >GaN^AlN, i.e. 
inversely with the cation size. While qualitatively simi- 
lar, our results differ somewhat from those of Ref. ^ for 
GaN. This discrepancy is probably due to the smaller 
cells and limited relaxations (first layer only) in Ref. 

In Table || the atomic displacements of first- and 
second-layer atoms are listed for GaN and GaAs 
(£=[001], z=[110]). For GaN the largest shifts are those 
of the surface cation, whereas the other atoms shift only 
slightly and almost rigidly. In GaAs, the displacements 
along z are much more relevant for both anion and cation. 
They move far away from each other, their vertical dis- 
tance [d=0.25 a (1 - 0.01 C B ) sin lo] being 0.69 A in 
GaAs against 0.23 A of GaN, on account of a more than 
double rotation angle, and of course of the 20 % larger 
lattice constant. 



III. (1010) WURTZITE SURFACES 

Wurtzitc is the most stable phase of III-V nitrides. Its 
(1010) surface is sketched in Fig. a. At variance with 
zincblende (110), symmetry only allows dimers rotation 
in the plane containing the [1010] and [0001] directions, 
i.e. ortogonally to the surface plane. Thus there is only 
one rotation angle: 9 — lo. The chemical picture closely 
resembles that of the (110). Instead of GaAs, we now 
consider ZnO (1010) as reference system. ZnO is one of 
the most ionic II- VI semiconductors, and it allegedly ex- 
hibits the same kind of relaxation anomalyCI under exami- 
nation here for the nitrides; thus, it represents a suitable, 
if extreme, case for comparison. 



Our results for the relaxations of wurtzite (1010) (Ta- 
ble |n|) basically confirm the findings for the zincblende 
(110) surfaces, with angles and bond contractions of the 
same order of magnitude (angles are somewhat smaller 
and, consistently, contractions are a bit larger). For GaN, 
previous calculationsu gave comparable, in fact somewhat 
smaller rotation angles. For ZnO our values can be com- 
pared with theoreticalcl and LEED resultsj!3 and are seen 
to agree well with the latter. Among the important fea- 
tures we note the close similarity between GaN and ZnO 
and the highly ionic character of A1N (see Sec. [v]); also, 
the "ionicity" trend InN^GaN^AlN is confirmed. 

In Table IV we list the atomic displacements in the first 



and second surface layer for GaN and ZnO. For GaN, the 
first-layer anions move upward, the cations downward. 
The separation along z is 0.36 A against 0.22 A of Ref. 

As a consequence our 8 is ~ 40% greater. Also, we 
find that even the second- layer cation moves upward size- 
ably, whereas in Ref. |§| changes in the second layer are 
moderate. Finally, for ZnO both surface atoms go down, 
but their distance along z (0.36 A) equals that of GaN. 



IV. SURFACE ENERGIES AND ELECTRONIC 
STATES 

In Table ^ we report the surface formation energies of 
all compounds studied, a is the surface energy per atom 
of the fully relaxed structure, Act the energy gained upon 
relaxation. Our results agree welLwith previous data 
for GaN (10T0)la and GaAs (110)£2 The formation en- 
ergy per atom may be roughly understood as the energy 
needed to break a single bond, i.e. 1/4 the cohesive en- 
ergy per atom: indeed, at least for the cases in Table [v|, 
a is close to i? co h/4. 

a is also reported in Fig. [| to make trends easily de- 
tectable. For the nitrides, the (110) surfaces energies are 
~60% larger than in GaAs. This difference is enhanced 
by relaxations, that strongly reduce the surface energy of 
GaAs. For the (1010) surfaces, energy differences range 
in an interval of ~ 0.2 eV. Noticeably, a grows along the 
same pattern of ionicity observed previously for bond 
contractions (InN^GaN^AlN). Finally, the (110) and 
(1010) surface formation energies of a given compound 
are similar, although the relaxation energy is larger for 
the latter. 

It is overall evident that as far as non-polar surfaces are 
concerned, the nitrides are closer to a highly ionic com- 
pound such as ZnO than to GaAs. Similar conclusions 
have been drawn from recent studies on spontaneous pa-, 
larization and piezoelectric constants of bulk nitrides 
indeed, all the data suggest that the nitrides are even 
more extreme in their deviation from typical III-V be- 
havior than ZnO compared to typical II- VPs. It is ap- 
propriate to check if such a behavior is also mirrored in 
the electronic properties. Indeed, while no surface states 
are wesent in the gap at the (110) surfaces of GaAs and 
GaPta because of dimcr rotation, for ZnO the occupied 
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dangling^bond surface state has been predicted to lay in 
the gap.u A previous calculations for GaN (1010) found 
the occupied anionic surface state to lay slightly (~0.1 
eV) below the valence band top. We find similar results 
(Fig. ^]), with the anionic surface state touching the va- 
lence top but still remaining within the band edge at 
r. Similar results are found for the other nitrides. In 
agreement with the detailed analysis of Ref. [9] for ZnO, 
the empty surface state corresponding to the remants 
of the cation dangling bond is prevailingly s-like. The 
filled surface state just above the valence band corre- 
sponds to the anion p-like back-bonds. It should also be 
mentioned that recent results by Hirsch et al£3 suggest 
that the occupied surface states just mentioned lay in- 
deed completely inside the gap. The difference to ours 
and other previous results should be attributed to an im- 
proved treatment of the semicore d states in Ref. |l8|. 



V. DISCUSSION 

To describe the relaxation mechanism, it is useful to 
consider separately three items. 

First, on all the (110) surfaces of binary tetrahedrally 
coordinated AjyB s _N compounds, a charge transfer oc- 
curs from the cation dangling bond into the anion dan- 
gling bond of the as-cleaved surface. This is a purely 
electronic-structure effect, occurring even at zero rota- 
tion angle: the cation dangling bond state is much higher 
in energy, and it fully transfers its electron occupancy 
into the anion state. 

Second, the surface dimer rehybridizes towards a cation 
sp 2 -like /anion p-like configuration. This entails a ro- 
tation of the dimer (a combined anion-upward, cation- 
downwards motion). This rotation is accompanied by a 
lowering of the energy of the occupied anion dangling 
bond state, and an increase of the energy of the cation 
empty state. This is precisely the reason why the cation- 
anion dangling bond occupation transfer is desirable for 
this rotation to happen. Since this rehybridization is 
qualitatively a kind of reverse of sp 3 hybridization, it is 
expected to be most favorable when the hybridization 
energy gain is low to begin with. 

Third, the charge within the surface dimer is asymmetric 
towards the anion, because of (a) chemical bond ionicity, 
in the spirit of e.g. the Garcia-CohenEJ charge asymme- 
try, and on top of that (b) the dangling bond transfer. 
Therefore the dimer rotation, with the ensuing anion dis- 
placement out of, and away from the cation plane, costs 
electrostatic energy. The energetic cost will be larger, 
the more asymmetric the charge distribution is (see e.g. 
Ref. ^). A key point is now that, in all materials, there 
is always a complete cation-anion dangling bond occupa- 
tion transfer: therefore, what matters is the net anion- 
cation charge-jasymmetry, that is largely equivalent to 
bulk ionicityEj The larger this is, the less the rotation is 
favored. To be precise, the dangling bond charge transfer 



will increase the local charge asymmetry (hence hinder 
rotation) more strongly in low-bulk-ionicity compounds: 
in the latter, the bulk (i.e. pre-cleavage) charge asym- 
metry is smaller than in more ionic compounds; thus, 
when the full dangling bond occupation is transferred 
to the anion, the net asymmetry increases more than in 
strongly ionic compounds. 

We can then rationalize the energetic balance as fol- 
lows. Rehybridization-plus-rotation is less costly when 
(a) the gain in the reverse process of sp 3 hybridization 
is low, i.e. qualitatively when the cohesive energy of the 
material is small, and (b) when the electrostatic cost of 
the outward rotation is low, i.e. when charge asymmetry 
is small, i.e. ionicity is low. The predicted trend is then 
that materials with small cohesive energies and ionicities 
will tend to have large rotations, and viceversa very ionic 
and strongly bound solids will tend towards small rota- 
tions. While "small" is to be understood in a relative 
sense, e.g. for GaAs compared to GaN, or ZnS com- 
pared to ZnO, the nitrides and ZnO are more ionic that 
all zincblende and wurtzite III-V's and II- VPs whatever 
the ionicity scale.E3E!J This picture agrees nicely with the 
calculated quantities for the nitrides as compared with 
other III-V's, asrwell as with those for ZnO as compared 
to other H-VI's;B both ZnO and the nitrides have both 
larger cohesion and ionicity (on any scale), and smaller 
rotation angles than their companion materials. (In the 
same direction, note that the dimer rotation can be inter- 
preted as a frozen-in zone-center A\ surface phonon;EU as 
all other phonons in the nitrides and ZnO, this mode is 
stiffer, hence more energetically costly, than in the other 
III-V's and II-VI's, respectively.) 

Indeed, large ionicities and cohesive energies, and 
hence small rotations are to be expected for first-row an- 
ions. The basic reason is that first-row atoms such as 
N and O have a very deep potential for the valence p 
states (and no p core-orthogonality constraint), whence 
stronger bonding and larger ionicity than with other 
cations. There is indeed a rather abrupt change in ro- 
tation angles (and in other properties too) for first-row 
anion both in II-VI's (see CdS vs ZnO in Ref. ||) and 
in III-V's; in this sense selenides are arsenidc-like, sul- 
phides are phosphide-like, and ZnO is GaN-like. Inter- 
estingly, a similar behavior is observed in the piezoelec- 
tric constants,Q which increase strongly as the anion de- 
creases in size. For the nitrides they are large and pos- 
itive, against the small and negative values of normal 
III-V's; for ZnO they are large and positive, against pos- 
itive and small in other II-VI's (the II-VI-III-V difference 
is due to a changed balance of the electronic and ionic 
components). 

Note however that one does not expect these trends to 
hold for any cation, in particular small ones. The trend 
for anions getting heavier and cations lighter is towards 
an effective exchange of roles (e.g. in boron compounds, 
this is reflected in anomalies of structural transitions un- 
der pressure), which blurs the picture somewhat. 

Dehybridization and charge asymmetry contrast each 
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other also in geometrical terms, i.e. the larger the layer 
rotation, the smaller the bond contraction. This can be 
seen by a simple geometrical argument. In Fig.|] an ide- 
alized picture of the surface profile is shown. The dashed 
line refers to the unrelaxed surface. If we keep the sur- 
face anion fixed and let the cation relax onto the plane 
formed by first- and second-layer anions (thus undergo- 
ing an ideal sp 3 — > sp 2 rehybridization), we have 9 ~35° 
and a bond contraction of ~5%. If the surface anion 
relaxes upwards, i.e. towards a more pure p-like config- 
uration (which indeed it does), then bond contraction 
tends to be suppressed. This is the case for GaAs, where 
large rotation angles (~ 30°) and small bond contractions 
(~ 1%) indicate that dehybridization dominates (also the 
case for other III-V's such as GaP). On the contrary, the 
small rotations in ZnO and the nitrides are accompanied 
by relatively large bond contractions, consistently with 
the more critical balance of electrostatic repulsion and 
deybridization. 

It is barely necessary to confirm explicitly that the 
nitrides are more ionic than other III-V's. Charge asym- 
metry increases with the electronegativity gap between 
cation and anion, commonly used as measure of com- 
pound ionicity. In Fig. ^ we report the experimental val- 
ues of electronegativity (i.e. one half the sum of atomic 
ionization potential and electron affinity) and the hard- 
ness (i.e. one half the difference of ionization and affinity) 
for the atoms under consideration. The main feature is 
that while cations behave quite similarly, this does not 
hold for anions, nitrogen having much larger values of 
X and r\. How does this influence charge asymmetry in 
compounds ? A semiquantitative estimate is-arovided by 
the electronegativity equalization principle,E3 which as- 
sumes the compound energy to be simply the sum of the 
atomic contributions. Upon compound formation, one 
obtains a charge transfer 



AN = 



Xb - XA 
2(r)B + tja) ' 



(1) 



which is depicted in Fig. ||. The charge transfer upon 
nitride formation is much lager, as a consequence of a 
greater ionicity. Use of other ionicity scalesf (Phillips, 
Pauling, etc.) will lead to the same qualitative concku. 
sions. For instance the charge asymmetry coefficientst£l 
g are 0.78, 0.79, and 0.85 for A1N, GaN, and InN re- 
spectively, a factor of ~ 2.5 larger than the 0.32 of GaN. 
Indeed, such huge difference is partially mitigated by the 
large chemical hardness of N in the denominator of Eq. 



VI. CONCLUSION 

In summary, non-polar surfaces of III-V nitrides pro- 
vide further evidence that the nitrides are closer to the 
extreme ionic limit (embodied e.g. by ZnO) than to nor- 
mal III-V compounds such as GaAs, in agreement with 



previous results on structural and polarization proper- 
ties. This strong ionic character causes the prevalence of 
dehybridization in determining surface relaxations to be 
less pronounced than in other III-V's. The same reason- 
ing applies to ZnO with regard to other II- VI compounds. 
In the final analysis, it is the nature of the nitrogen anion, 
in particular its strong valence potential and the absence 
of core p states, that sets the nitrides apart from the other 
III-V's, just as the analogous properties of oxygen cause 
the major differences of ZnO and other II- V compounds. 



ACKNOWLEDGMENTS 

The calculations were performed on the IBM SP2 of 
CRS4 (Centro Ricerche, Sviluppo, e Studi Superiori in 
Sardegna). V.F.'s stay at WSI was supported by the 
Alexander von Humboldt-Stiftung. V.F. thanks Paolo 
Ruggerone for enlightening discussions and for reading 
the manuscript. 



1 F. A. Ponce and D. P. Bour, Nature 386, 351 (1997). 

2 Diamond, Silicon Carbide and Related Wide Ban 
Semiconductors, edited by J. T. Glass, R. Messier, and N. 
Fujimori, MRS Symposia Proceedings No. 162 ( Materials 
Research Society, Pittsburgh, Pa, 1990); M. J. Paisley, Z. 
Sitar, J.B. Posthill, and R. F. Davis, J. Vac. Sci. Technol. 
A. 7, 701 (1989); G. Martin, S. Strite, J. Thornton, and H. 
Morkoc, Appl. Phys. Lett. 58, 2375 (1991). 

3 P. Manca, J. Phys. Chem. Sol. 20, 268 (1961). 

4 V. Fiorentini, M. Methfessel, and M. Schemer, Phys. Rev. 
B 47, 13353 (1993). 

5 F. Bernardini, V. Fiorentini, and D. Vanderbilt, Phys. Rev. 
Lett. 79, 3958 (1997). 

6 F. Bernardini, V. Fiorentini, and D. Vanderbilt, Phys. Rev. 
B 56, R10024 (1997). 

7 J.E. Jaffe, R. Pandey, and P. Zapol, Phys. Rev. B 53, 4209 
(1996). 

8 J. E. Northrup and J. Neugebauer, Phys. Rev. B 53, 10477 
(1996) 

9 P. Schroer, P. Kriiger, and J. Pollman, Phys. Rev. B 49, 
17092 (1994); J. Pollman, private communication. 

10 J. E. Jaffe, N. M. Harrison, and A. C. Hess, Phys. Rev. B 
49, 11153 (1994). 

11 R. Dreizler and E. K. U. Gross, Density- functional theory 
(Springer, New York 1990). 

12 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 556 
(1980). 

13 J. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

14 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

15 See e.g. G.-X. Qian, R. M. Martin, and D. J. Chadi, Phys. 
Rev. B 38, 7649 (1989). 

16 F. Bechstedt and R. Enderlein, Semiconductor Surfaces 
and Interfaces (Akademie-Verlag, Berlin 1988) 



4 



A. Paxton, and P. Mark, Phys. TABLE [II Dimer rotation angle (9) and relative bond 

1010) surfaces. Labels a and 6 refer 
| and ^, respectively. Experiments 



17 C. B. Duke, R. J. Meyer ; 
Rev. B 18, 4225 (1978). 

18 G. Hirsch, P. Kriiger, and J. Pollmann, unpublished. 

19 A. Garcia and M. L. Cohen, Phys Rev. B 47, 4215 (1993); 
ibid., 4221. 

20 W. Monch, Semiconductor Surfaces and Interfaces 
(Springer Series in Surface Science 26, Springer, Berlin 
1993). 

21 Y.R. Wang and C.B. Duke, Surf. Sci. 205, L755 (1988); P. 
Santini, L. Miglio, G. Benedek, U. Harten, P. Ruggerone, 
and J. P. Toennies, Phys. Rev. B 42, 11942 (1990). 

22 R. G. Parr and W. Yang, Density- Functional Theory of 
Atoms and Molecules (Oxford University Press, New York, 
1989). 



contraction (Cb), for ' 
to results from Rcfs. 
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TABLE | Surface dimer rotation angles 9 and u> (see Fig. 
||), and relative bond contraction Cb for zincblende (110) 
surfaces. 9 a , uj a and Cg are from Ref. ft]. For GaAs ex- 
perimental values are also shown (Ref. H3). 





A1N 


InN 


GaN 


GaAs 


e 


11.7° 


14.4° 


14.3° 


30.1° 








2.06° 


24.3° 


^Expt. 








31.1° 




5.8° 


7.4° 


7.3° 


16.5° 


w a 






1.0° 


13.4° 


^Expt 








16.7° 


C B 


5.9% 


4.3% 


4.9% 


0.9% 


C% 






6.5% 


1.3% 


QExpt 








2% 



TABLE H Displacements from ideal positions (in A) of 
anions {An) and cations {Cn) in first and second layer of 
GaN and GaAs (110). £=[001] and £=[110]. 





GaN 




GaAs 






Ax 


Az 


Ax 


Az 


Am 


-0.04 


0.05 


0.15 


0.42 


Cm 


0.17 


-0.18 


0.37 


-0.27 


An 2 


-0.05 


0.02 


0.06 


0.13 


Cn 2 


-0.03 


0.07 


0.08 


0.23 



TABLE [V Atomic displacements of first (An\ and Gn\) 
and second {An?, and CV12) layer from ideal positions (in 
A) for theJIOlO) surface of GaN and ZnO. (£=[0001] 
and £=[1010]). An and Cn indicate anion and cation, 
respectively. Superscript a indicates results from Ref. |^. 







GaN 






ZnO 






Ax 


Ax* 


Az 


Az* 


Ax 


Az 


Am 


0.04 


0.01 


0.08 


0.02 


0.02 


-0.13 


Cm 


-0.15 


-0.11 


-0.28 


-0.20 


-0.14 


-0.50 


An 2 


0.04 


0.05 


-0.02 


0.05 


-0.02 


-0.09 


Cn 2 


0.05 


0.05 


0.15 


0.05 


0.03 


-0.09 



TABLE ^ Surface formation energies (ct), relaxation en- 
ergies (Act) and cohesive energies per bond (i.e. £' co h/4, 
where -E" co h is the cohesive energy per atom) . Results are 
in eV/atom. 
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FIG. 1. Side view of (110) ZB surface for III-V semiconduc- 
tors. White spheres are anions (An), black ones are cations 
(Cn). 



FIG. 4. Formation energies a of (110) and (1010) surfaces. 




FIG. 2. Dimer rotation at (110) surface; 9 and uj are two 
independent parameters. 
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FIG. 5. Side view of the (110) surface. Dashed line de- 
notes the ideal structure, full lines two possible atomic rear- 
rangements, one with the anion kept fixed in its ideal position, 
the other with anion shifted upward. 
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FIG. 3. Side view of the relaxed (1010) surface, 
spheres are anions (An), black ones cations (Cn). 
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FIG. 6. Band energies of GaN (1010), plotted for k run- 
ning along the diagonal of IBZ (shown in the inset). Black 
circles are anionic (Sn) and cationic (Sgo) surface states. 
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FIG. 7. Electronegativity x and hardness r\ for the atomic 
components of some III-V semiconductors. N is considerably 
more electronegative than the other anions. 



FIG. 8. Charge transfer per dimer given by the electroneg- 
ativity equalization model. 
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